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Abstract 

Computations of the Julia and Mandelbrot sets of the Riemann zeta func- 
tion and observations of their properties are made. In the appendix section, a 
corollary of Voronin's theorem is derived and a scale-invariant equation for the 
bounds in Goldbach conjecture is conjectured. 

1 Introduction and Motivation 

In this paper, we shall present a series of diagrams showing the Julia set of the 
Riemann zeta function and its related Mandelbrot set. This computation of the 
Julia and Mandelbrot sets of the Riemann zeta function may be, with reference 
to the existing literature, the first attempt ever made. 

As we shall see, the Julia and Mandelbrot sets of the Riemann zeta function 
have unique features and are quite unlike those of any elementary functions. 

We shall also point out, in the Appendix, two known observations of ap- 
proximate self-similarity in Number Theory — in the images of Riemann zeta 
function itself, and in the problem of the Goldbach conjecture where a func- 
tional equation or scale-invariant equation for the ds in Goldbach conjecture is 
conjectured. 



2 Properties of the Iterated Maps of the Riemann 
Zeta Function 

Definition 1 (Iterated Maps of the Riemann zeta function) 

Denote the iterated maps of the Riemann zeta function £(s) as Cn '■ C i— > C 
such that 

Cn-.s ^ c(C(---COO---)) (i) 

V v ' 

n-times 



Definition 2 (Attractor and Repellor) 

A fixed point p of f(x) is an attractor if \ f'(p)\ < 1, a repellor if |/'(p)| > 1. 
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Lemma 1 (An Asymptotic Property of C(s)) 

log|C(-|a|)| = Hlog(|a| + l)-H(l+log27r) + 0(log|(7|) (a G R) . (2) 
Proof 

From the functional equation of the Riemann zeta function 

C(l- S ) = 2(27r)-*r( S ) cos(^) C(-), (3) 

we have 

\((a)\ < |(2vr)^- 1 r(l-a)C(l- ( r)| 

as | sin(o"7r/2)| < 1. 

As a — > oo, C(°") ~~ *■ 1 faster than (1 + o~~ q ) — > 1 for any positive real 

C(i+|(x|) = i + o(H-«), 
io g (c(i+M) = o(H-«), 

iog(r(i+|(r|)) = (i + H)iog(i + H)-(i + H) + o(io g H) 

(by Stirling's formula) . 

Thus, 

log|C(-k|)| = -Hlog(27r) + (l + H)log(l + H)-(l + H) + 0(logH) 
= |CT|log(H + l)-|<7|(l + log27r) + 0(logH) . 

□ 

Theorem 1 (Fixed Points of £ n ) 

The Riemann zeta function £(s) has one attractor fixed point at 

a = a = -0.2959050055752 . . . 
but infinitely many repellor fixed points. 

Proof 

For a > 0, y = a intersects y = £(er) at only one point, a = a + = 
1.833772 .... C(«+) is thus a fixed point of ( n . Since |C'(a+)| = 1-3742 . . . > 1, 
the fixed point at a = q+ is a repellor. 

For -1 < a < 0, y = a intersects y = C(er) at a = a = -0.2959050055752 . . . . 
C(ot) is thus a fixed point of Q n . Since |C'(aOI — 0.51273 . . . < 1, the fixed point 
at a = a is an attractor. 

For a < —1, £(2(7 + 1) = at a G Z~, and C(°") oscillates with amplitude 
3> |<r| (by Lemma |l]). Thus for a < —1, y = o~ intersects y = C(°~) twice for the 
interval of every two zeros of C(°0> an d |C'(°") I ^ 1 f° r °~ m the neighbourhood of 
these intersections. C, n has thus infinitely many repellor fixed points for a < — 1. 
□ 
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3 Julia set and Mandelbrot set of ((s) +a for s,a G C 



For the quadratic map F(a, s) = s 2 + a where s, a £ C, we have the following 
definitions. 

Definition 3 (Julia set) 

The Julia set J of a map F is the boundary of the domains of the attractors of 
the iteration of F. 

For example, if A\, A 2 are attractors of the iteration of a map F, D(A\), D(A2) 
are the domains of the attractors respectively, and dD(A\), dD(A2) are the 
boundaries of the domains respectively, then the Julia set of F is J = dD{A\) = 
dD(A 2 ). 

Definition 4 (Mandelbrot set) 

The Mandelbrot set M of a map F(a,s), for a chosen initial iteration value 
s = sq, is the set of values of the complex parameter a which, in the iteration of 
F(a, s) with the chosen initial value sq, does not lie in the domain of attraction 
of complex infinity. 

M = {a : a G C, lim F n (a,so) 00} . (4) 
The Mandelbrot set is related to the Julia set in the following theorem. 

Theorem 2 (The Julia- Fatou Theorem) JJ, p. 105] 
The Julia set J of a map F(a, s) for the given value of a is connected if and 
only if lim F n (a, 0) -A 00, i.e. if and only if a belongs to the Mandelbrot set M 

of the map F(a, s) for s = 0. 

Though the above definitions are for polynomials and are not the most 
general for non-polynomials, we shall choose to adopt these definitions here to 
study, for starters, a particular one-parameter family of the iterated maps of 
the Riemann zeta function, i.e., that of F(a, s) = £(s) + a. 

The Riemann zeta function £(s) has one attractor at fixed point s = a (by 
Theorem ||). It is not known if Q{s) has any non-real fixed point, and this 
remains to be proved. 

The Julia set of £(s) has been numerically computed with Mathematica with 
the algorithm downloadable from the Internet at 

http : / / www, damtp. cam .ac.uk /user /scw2 1 / riemann / 

The results are presented in Figures [l], ^, ^ and This may be the very 
first computation of the Julia set of £(s) ever attempted. 
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Re(a) 

Figure 1: The Julia set of ((s) for -20 < Re(s) < 10 and < lm(s) < 30 . 



The black region represents the basins of attraction of the attractor of the 
fixed point at s = a and those of attracting cycles. The other attractor is at 
oo. The Julia set is the boundary of the black region. 

The regions of different shades represent the various rates of amplitudes 
in iteration — > oo, with the lighter shades having slower rates, and the darker 
shades, faster rates. 

There is a reflection symmetry about the Re(s)-axis due to the complex 
conjugate property of an analytic complex function. 



-20 -15 -10 -5 5 10 




-20 -15 -10 -5 5 10 



Re(a) 

Figure 2: The Julia set of ((s) for -20 < Re(s) < 10 and 30 < \m(s) < 60 . 

Although a "piece" of the black region up along the lm(s)-axis may resemble 
some other nearby "piece", when examined in detail, each "piece" is unique and 
is distinct from any other "piece" . 
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-20 -15 -10 -5 5 10 




-20 -15 -10 -5 5 10 



Re(a) 

Figure 3: The Julia set of ((s) for -20 < Re(s) < 10 and 60 < \m(s) < 90 . 

The "pieces" of the black region seem to run indefinitely up the lm(s)-axis 
and lies on Re(s) > thereon. 
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-8.5 -8.25 



-8 -7.75 

Re(a) 



-7.5 -7. 25 



-7 



Figure 4: Zooming in to reveal the self-similarity of the Julia set of £(s) for 
-8.5 < Re(s) < -7 and 6.5 < lm(s) < 8 . 

We shall now turn to the Julia sets of F(a, s) = ((s) + a for a € C. We shall 
evolve the Julia sets of F{a, s) from one value of a to another along a path in 
the complex plane of a, and study the dynamics in the change of the structure 
of the sets. In particular, we shall evolve the Julia sets of F(a, s) along a = 
to a = — 1 along the real axis of a as in Figure [|, and similarly but in opposite 
direction from a = to a = 1 along the real axis of a as in Figure [6|. 

A numerical observation is that if sequence F n (a,s) in the Julia set itera- 
tions converges, then F n (a,s) tends to a fixed point which is a function pj(a) 
independent of s. This translates to the following conjecture: 

Conjecture 1 

Given a fixed a, + a = s has either one or no finite solution. 
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(iv) a = -0.1 (v) a = -1 



Figure 5: Julia sets of F(a, s) = ((s) + a 
for -20 < Re(s) < 10 , < lm(s) < 30 , and 
(i) a = 0, (ii) a = -0.001, (iii) a = -0.01, (iv) a = -0.1, (v) a = -1. 
The horizontal axis is Re(s), and the vertical axis is lm(s). 
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(iv) a = 0.1 (v) a = 1 

Figure 6: Julia sets of F(a, s) = ((s) + a 
for -20 < Re(s) < 10 , < lm(s) < 30, and 
(i) a = 0, (ii) a = 0.001, (iii) a = 0.01, (iv) a = 0.1, (v) a = 1. 
The horizontal axis is Re(s), and the vertical axis is lm(s). 
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We shall now consider the Mandelbrot set of the Riemann zeta function. 
Given the map F(a, s) = £(s) + a for which the initial iteration value s = sq is 
a zero of the Riemann zeta function C( s )) a value of the complex parameter a 
belongs to the Mandelbrot set of F if the iterated value F n (a, so) does not tend 
to complex infinity in the limit the order of iteration, n — > oo. 

The results of the computation of the Mandelbrot set of F(a, s) = £(s) + a 
where the initial iteration value s = sq E {s : ((s) = 0}, are presented in Figures 
[7|, H and ^. This may also be the very first attempt at such computation. 

However, the following about the Mandelbrot set of the Riemann zeta func- 
tion was observed numerically and remains to be proved: 

Conjecture 2 

For the map F(a, s) = £(s)+a, if a belongs to the Mandelbrot set, then F n (a, so) 
in the Mandelbrot iterations tends to a fixed point which is a function pm{o) 
independent of sq, or tends to a n-cycle where n G Z + . 

In addition, if Theorem ^ which applies to the quadratic map also applies 
to F(a,s) = C(s) + a, from Figure [7], we can see that a = belongs to the 
Mandelbrot set of F(a, s) = ((s) + a for s = 0, we deduce that the Julia set of 
((s) is connected. 
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-28 -23 -18 -13 -8 -3 2 



Re(s) 

Figure 7: Mandelbrot set of F(a, s) = ((s) + a, 
for which s is a zero of C( s )> 
-28 < Re(a) < 2 and < lm(o) < 30 . 



The black region represents the Mandelbrot set. The regions of different 
shades represent the various rates of amplitudes in iteration — » oo, with the 
lighter shades having slower rates, and the darker shades, faster rates. There 
is a reflection symmetry about the Re(s)-axis due to the complex conjugate 
property of an analytic complex function. 

The overall picture has the feature of a "comb" or "bristles of a brush" with 
"hairy ends". It is interesting to note that the "bristle" structures resemble 
those of the Julia sets for \a\ » in Figures |6] and [5] although the structures 
at the end of the "bristles" are different. 
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-5.5 -4.5 -3.5 -2.5 -1.5 -0.5 0.5 



Re(a) 

Figure 8: Zooming in to reveal the self-similarity 
of the Mandelbrot set of F(a, s) = ((s) + a 

for which s is a zero of C( s )> 
-5.5 < Re(a) < 0.5 and 18<lm(a)<24. 

The picture zooms onto the "hairy end" of a "bristle" only to reveal more 
"bristles" structures with "hairy ends" features, and ad infinitum. 
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Re(a) 

Figure 9: Mandelbrot set of F(a, s) = £(s) + a for s = 0, 
-28 < Re(a) < 2 and < lm(a) < 30 . 

We can see that a = lies in the black region, and thus a = belongs to 
Mandelbrot set of F(a, s) = ((s) + a for s = 0. 

We shall consider the more general case, the iterations of fully-parametrized 
family of maps of the Riemann zeta function, in the next paper. 
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Appendix 

A Examples of Approximate Self-Similarity and Scal- 
ing in Number Theory 

A.l Approximate self-similarity in the Riemann zeta function 

We shall define approximate self-similarity and show that the following theorem 
due to Voronin implies that there are approximate self-similarity in the images 
of the Riemann zeta function. 

Definition 5 (Approximate Self-Similarity to order (e, r)) An n- dimensional 
map f has approximate self- similarity to order (e, r) when given e > and r > 0, 
there exist X a ,Xi, and r& such that \f(X — X a ) — f{ri{X — Xf,))\ < e for all 
\X\ < r, where X = (x\, X2, • • • , x n ). 

Theorem 3 (Voronin) J]]/ 

Let < r < 1/4 and let f(s) be a complex function analytic and continuous for 
\s\ < r. If f(s) ^ ; then for every e > 0, there exists a real number T = T(e, /) 
such that 

max |/(s) - C(s + 3/4 + %T) \ < e . (5) 

\s\<r 

Theorem 4 There exists approximate self- similarity to order (e, r) for any 
given e > 0, < r < 1/4 in the images of the Riemann zeta function. 

Proof 

Choose f(s) = C,(pe t9 s + a) with fixed p, 6 € R and fixed a £ C where 
p > 0, 9 G [0, 2tt), < |s| < 1/4, and choose e to be arbitrarily small. The 
expression pe %e s + a with fixed p,9,a describes a disc centred at a with radius 
p/4 and orientation rotated at angle 9. 

By the inequality || in Theorem we have self-similarities in the Riemann 
zeta function up to the order of e between the images of the discs at differ- 
ent centres, scales and orientations (recentred by choosing a different fixed a, 
rescaled by choosing a different fixed p, and rotated by choosing a different 
fixed 9), and the images of the discs along the upper half of the critical strip at 
s + 3/4 + iT with different values of T. □ 
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A. 2 Scale-invariant equation for bounds in Goldbach conjec- 
ture 

Conjecture 3 (Goldbach conjecture) 

Every even integer 2n > 4 can be written as a sum of two primes. 



Definition 6 

Define G(2n) as the number of distinct decompositions of the even integer In 
into a sum of two primes (3+5 and 5+3 are counted as identical decomposition). 

Define Gl{x) and Gjj(x) as the lower and upper monotonic bounding curves 
respectively of G{2n) such that 

dGjJjK) dGujx) 

: > and : for all x , 

dx dx 

and Gr.(2n) < G(2n) < Gjj(2n) for all even integer 2n > 4. 

Thus, Goldbach conjecture states that Gl(2u) > 1 for all even integer 
2n > 4. However, the plots of G(2n) in Figure |l^ seem to show that the shapes 
of the bounding curves are invariant under scaling. This observation leads us 
to suggest a new conjecture: 



G{n) 




n n n 

Figure 10: Scale invariance in the bounding curves of G(2n) and the approxi- 
mate self-similar pattern in G{2n) 



Conjecture 4 (Scale- invariant equation for bounds in Goldbach conjecture) 

Gl{x) satisfies the scale invariance equation 

G L (ax) G L {bx) 
G L (a) 



G L {b) 

and similarly for Gu(x). 



G L {bx/a) = G L {x) , (6) 



Although f{ax) = f(a)f(x) (for all a) satisfies (||), the form of Gl{x) is 
not in the form of f{x). Since the simple power law y(x) = (3 x x , (3, A E R is in 
the form of f(x), we thus deduce that Gl{x) is not a simple power law. 
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The patterns in G(2n) as illustrated by Figures O and 12 are due to the dis- 
tribution of prime numbers. Hardy and Littlewood ||] had derived an estimate 
for G{2n) which accounts for the patterns, 

G(2n) ~ -77 

V ' logn 2 { L p-2 

y ° ' p>2, p\n^ 

where p runs over all prime numbers, p\n denotes that the product is taken only 
if n is divisible by p, and C is a constant, 

c -n( i -^n ? )- M « w -- 




Figure 11: The subset of {G(2n)} with n being a multiple of 3, {G(2(3n))}, is 
labelled by blue dots. The subset which is not, {G(2(3n + l))}n{G(2(3n + 2))}, 
is labelled by red dots. 



G(n) 




n 



Figure 12: {G(2n) \ n = mod 3, n = mod 5} is labelled by magenta dots, 
{G(2n) | n = mod 3, n ^ mod 5} by blue dots, {G(2n) \ n ^ mod 3, n = 
mod 5} by green dots, {G(2n) \ n ^ mod 3, re ^ mod 5} by red dots. 



However, the scale- invariant property of the bounding curves Gl(x) and 
Gjj(x) of G{2n) as expressed in (||) is a new observation. 
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A. 3 Open Problem 

1. Derive the exact expressions for the bounding curves Gl{x) and Gu{x). 

We shall attempt to derive Gl{x) and Gjj(x) in the next paper. 
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